Brain Topography

Background: Electroencephalographic (EEG) potentials evoked by transcranial magnetic stimulation (TMS) offer a direct window into cortical dynamics. Yet, a systematic exploration of their morphological features, analogous to sensory-evoked potentials, is lacking, especially for stimulation outside the motor cortex. Aim: To obtain region-specific properties of frontal, parietal and occipital networks from the time course of TMS-evoked potentials (TEPs). Materials and Methods: We implemented and applied an automatic procedure to compute peak-to-peak amplitude, peak latency, and inter-peak interval of TEPs recorded from 40 neurotypical subjects stimulated over left occipital (n=25), parietal (n=25), and frontal (n=25) cortices. Results: Occipital TEPs showed the largest peak-to-peak amplitude and longest latency of the first waveform component, independently of stimulation intensity and consistent with the recruitment of a large patch of densely interconnected neurons. Concerning later components, both latency and inter-peak interval systematically decreased along the posterior-to-anterior axis, reflecting progressively faster recurrent dynamics from the alpha-dominated occipital circuitry to the tightly coupled loops between frontal cortex and subcortical structures. Parietal TEPs showed intermediate amplitude and latency measures, consistent with the heterogeneous cytoarchitectonic and connectional organization of the superior parietal cortex. Conclusions: Our findings suggest that TEP morphology is shaped by the distinct properties of the stimulated networks, with early amplitude reflecting the extent of local recruitment and later temporal features tracking the rhythm of recurrent activity. This work offers a mechanistically grounded and practically accessible approach, also released as a Python-based tool, that allows to characterize cortical reactivity across different brain-states and populations.

Learning new information in the presence of distracters and changing conditions requires the ability to adapt. In the brain, this adaptive capability has been linked to dynamic interactions between attention and working memory, which enable the selective filtering of irrelevant input while preserving behaviorally relevant information. Specific neural oscillations have been implicated in this process. Here, we introduce a phenomenological data-driven framework for oscillatory network modeling that learns condition-dependent coupling laws directly from neural recordings and enables inference of condition-dependent directed pathways. We apply our approach to magnetoen-cephalography (MEG) data collected while participants performed a working-memory task with and without distracters. Recall dynamics in the non-distracter condition are first modeled using a linear oscillatory network in which each region of interest is represented by two alpha-band harmonic oscillators. We use universal differential equations (UDE), an extension of neural differential equations, to capture distracter-induced changes in coupling laws. Symbolic regression is then used to interpret the modifications identified by UDE as nonlinear functions, and an additional method is proposed to identify the directed pathway from the newly emerging nonlinear terms in the dynamics of brain regions of interest. Despite inter-subject variability, working memory recall data from all four participants examined under distraction showed the emergence of a pathway from the dorsolateral prefrontal cortex (dlPFC) to the primary visual cortex (V1). This finding is consistent with the established role of the dlPFC in cognitive control and suggests that distracter processing recruits a directed interaction from prefrontal to visual regions. More broadly, our results illustrate that combining linear models whose parameters are learned from the data with universal differential equations augmented by interpretability methods enables the identification of condition-dependent coupling laws, their representation as interpretable mathematical functions, and the discovery of candidate directed pathways underlying adaptive changes in oscillatory networks without requiring strong prior assumptions about the underlying mechanisms.

Derived from previous observations on equal and cross-frequency coupling, we evaluated the proposal that equal and cross-frequency phase synchronization may characterize the integration-segregation perspective of cerebral sensory-motor processing. Using brain recordings obtained in normal conditions and in conditions of diminished sensory input (eyes closed wakefulness, sleep and coma, when there is presumably less functional segregation of sensory-motor processing in neural networks), we assessed potential differences in partitioning of the synchrony state space linked to cross-frequency synchronization. More partitions were found in conditions of decreased sensory input. In addition, there was a less complex synchrony state space in cross-frequency as compared with equal-frequency coupling, in terms of fewer connectivity configurations. These results support the idea that equal-frequency coupling favours integration from multiple brain regions occurring in a complex synchrony state space rich in possible connectivity configurations, whereas cross-frequency coupling contributes to segregation, or localized sensory-motor transformations taking place in specific brain areas. This evidence may contribute to new considerations about the much-discussed role of multi-frequency relations in neuronal activity, and how the structural and functional modular organization of the nervous system is able to generate the coordinated activity needed for conscious and appropriate cognitive behaviors in complex environments.

Transcranial Electric stimulation (tES) is a safe and noninvasive technique increasingly used in treating brain disorders. Despite many studies on tES, there is still a lack of understanding of its mechanisms at the microscale network level. This is crucial for optimized parameter selection in therapy approaches such as the treatment of pharmacoresistant epilepsy. In this study, we made use of a recently published neuroinspired microscale model of the neocortex, known as NeoCoMM, and integrated a "Lambda E"-based model of tES. This updated version was used to investigate the acute effects of tES (tDCS and tACS), on the neural activity of various neuron types in both healthy and epileptic brain states. Results showed that in the case of healthy alpha and gamma rhythms, tACS induced electric field entrainment at the peak power frequency of the network oscillations as measured by the Local field Potentials (LFPs). This resonance-like entrainment was independent from the individual firing rate of cell types. For epileptic activity, tACS did not provide consistent results. Cathodal tDCS resulted in a promising decrease in hyperexcitable activity throughout simulations. These results advance our understanding of the impact of tES on network dynamics at both the extracellular and intracellular activity levels. Author summary

The behavioural activity of freely moving animals is a confounding factor that affects the recording, analysis, and final results of animal EEG experiments. Along with the lack of standardisation in animal in vivo electrophysiology experiments, this could lead to huge inconsistencies, especially in the analysis of centrally acting drugs. Therefore, the main aim of this paper is to investigate the effects of behavioural activity versus inactivity on the multichannel EEG in freely moving rats. In a large sample (n = 116) of waking recordings from 12 cortical electrodes (ECoG) in Wistar rats, we evaluated behavioural activity-related changes in the power spectrum, current source density, and power-based global functional connectivity (GFC) in a 3D rat brain model, according to the TOHOKU Rat Brain Atlas. The main findings were that behavioural activity induced 1) a robust power increase in 6-8 Hz, peaking at 7 Hz with maximum changes over the parietal and temporal cortex, 2) an increase in gamma power (30-80 Hz) across the whole brain, 3) a decrease in delta (1-4 Hz) and beta (12-30 Hz) power across the whole cortex. Changes were also localised in subcortical regions, particularly in the diencephalon/thalamus. The GFC analysis showed a similar pattern of power changes across the 6-8 Hz, delta, and beta bands; however, GFC in the gamma band decreased. Again, the GFC analysis revealed changes in connectivity within subcortical structures, primarily in the thalamus. None of the measures was affected in the alpha band (8-12 Hz). These findings emphasise behavioural state as a critical factor influencing EEG outcomes, with important implications for the standardisation and translational validity of preclinical neurophysiological studies.

Intracranial neurophysiology studies have typically ignored signals from electrodes located in white matter (WM), assuming that their information content is artifactual or related to nearby gray matter (GM). Here, we tested the electrophysiological and functional features of signals recorded from different WM locations. Signals were recorded from 19 patients undergoing intracranial monitoring for drug-resistant epilepsy by means of stereo-electroencephalography (sEEG). Each sEEG electrode was classified into WM or GM based on the surrounding tissue. We obtained recordings from a total of 1,717 sEEG electrode contacts, 36% in WM, while the patients were in awake resting state (5 minutes). For each sEEG electrode, we employed a model-based spectral decomposition to separate periodic and aperiodic components, and we computed signal complexity metrics. For a subset of participants, we computed WM structural information from diffusion-weighted magnetic resonance imaging and we evaluated functional signals during a cognitive control task. Our results show that signals recorded from WM have different spectral features and higher complexity than GM. Complexity correlates positively with fractional anisotropy, and modulations related to behavior during the task were detected in WM. Overall, this indicates that WM signals carry information that may reflect signal propagation across WM fiber tracts.

Intense aesthetic experiences are among the most complex responses arising from the interaction of mind, brain, and context. Observations from fMRI suggest that when viewers feel highly moved by artworks, the underlying neural states differ from those accompanying less intense responses, particularly through recruitment of the DMN. Using electroencephalography and Bayesian category-specific cumulative link mixed models, we investigated whether such putative peak aesthetic responses exhibit threshold-specific neurodynamics rather than linear scaling with intensity. Twenty-two Chilean participants viewed 113 diverse local artworks whilst rating how moved they felt on a four-point scale. We analysed both canonical oscillatory power ({theta}-{gamma}) and aperiodic components (offset and exponent) during the contemplation window and the post-elicitor window. Threshold-specific effects were found: spectral features differentiated the highest rating category from moderate responses, rather than scaling uniformly across all intensity levels. During artwork visualisation, power in the {beta}1 and {beta}2 bands, as well as the interaction of {beta}1 with the 1/f exponent, predicted the transition to the most intense response; during the post-elicitor window, the aperiodic 1/f exponent predicted the transition from very low to higher-intensity responses. Modelling individual differences in spectral signatures (in the and {gamma} bands) credibly improved predictive performance (approximate leave-one-out cross-validation; elpd_loo), suggesting that neural variability reflects meaningful mechanistic heterogeneity in aesthetic processing rather than mere noise. These findings speak to a broader question, how the brain marks the intensity of conscious experience, and, more specifically, support the hypothesis that being intensely moved constitutes a qualitatively distinct neural state, characterised by specific configurations of oscillatory dynamics and cortical excitability that modulate the transition from low and moderate to peak engagement.

Substantial progress has been made in recent years on understanding how the human brain parses and processes natural speech. Much of this progress has been based on modeling how brain activity relates to the different acoustic and linguistic features of speech. By fitting and testing models based on those features, one can test hypotheses about the kinds of computations and representations the brain uses to convert speech sounds into understanding. While much of this work has focused on modeling BOLD activity using functional neuroimaging or intracranially recorded electrophysiological signals, the approach has also proven useful with MEG and EEG. Indeed, noninvasive EEG has certain advantages for studying speech processing in terms of translational research and application. Research over the last decade or so has shown that EEG can be successfully modeled based on numerous acoustic, linguistic, and paralinguistic speech features. However, an important unanswered question hangs over all of this work: namely, what constitutes a good model of EEG responses to natural speech? Or, to put it another way, how much variance in EEG recorded during natural speech listening is explainable as having derived from that speech input? The present study aims to tackle this issue. We do so under the assumption that the best model for a person's EEG response to natural speech is a set of EEG responses from other people listening to the same speech. Using this assumption, we construct inter-subject models using EEG from 19 healthy adult native speakers of English who all listened to the same audiobook. The model for each subject involves predicting their EEG data using (dimensionality-reduced) EEG from different numbers of other subjects and then extrapolating to estimate the total explainable variance in the target individual's response to speech. Following this, we show that linear models (temporal response functions) based on several commonly used acoustic and linguistic speech features can predict most - but importantly not all - of the estimated total explainable variance in EEG responses across subjects.

The brain requires coordination among different regions to execute cognitive tasks, which may involve both positive- and negative-correlations. The topology of these correlations may indicate the mechanism underlying brain functioning in a given state. Here, we study changes in the functional connectomes (FCs) of both the positive and negative-correlations across various cognitive task states relative to the resting state, using publicly available electroencephalographic (EEG) data. Considering the EEG-specific topographical cortical regions as topographical modules (TMs), we find that the FC comprising positive correlations (G+) is modular. In contrast, networks of negative-correlations (G-) are anti-modular, with more connections between TMs than within them, and are associated with improved overall topological efficiency. These functional networks also show variability across frequency bands and brain states. In the low-frequency delta band, resting states exhibit higher modularity and anti-modularity than task states; in contrast, in the high-frequency Gamma band, modularity and anti-modularity are much higher during task states than in the resting state. The k-core analysis of all networks further reveals differences: G+ is more hierarchical and robust than G- across all states. Moreover, the task-state networks are always more hierarchical than the resting-state networks across all frequency bands. In the high-frequency gamma band, they are also significantly more robust than the resting-state networks. These networks also differ in the topology of their innermost core constituents: the innermost core regions of G+ are randomly connected and spatially localized, mostly in posterior brain regions across subjects, in the high-frequency gamma band. Whereas those in G- are spatially de-localized, cover the extreme anterior and extreme posterior brain regions, and remain anti-modular in all the frequency bands. Overall, our analysis reveals the presence of an anti-modular organization of functionally specialized TMs alongside their modular organization and points to task- and resting-state differences in their topologies.

Music enjoyment and familiarity are closely related but often confounded in studies of neural music processing. Here, we investigated their distinct contributions to cortical oscillatory activity and neural tracking of music using electroencephalography (EEG). Thirty-two participants listened to self-selected all-time favourite songs, recent favourite songs, and tempo-matched songs from disliked genres. This novel paradigm dissociated familiarity from enjoyment by including highly enjoyed songs that differed in familiarity. Spectral power and cortical tracking (using Mutual Information) were analysed using linear mixed-effects models with enjoyment and familiarity ratings. Familiarity was associated with increased left-frontal alpha power, whereas enjoyment predicted increased theta and beta power, demonstrating distinct oscillatory signatures for these dimensions. An interaction revealed that the positive relationship between enjoyment and theta power was strongest for highly familiar music. Cortical tracking analyses showed that greater enjoyment was associated with reduced delta-band tracking, with a significant interaction indicating that this negative relationship was present for highly familiar songs but not for less familiar songs. These findings indicate that enjoyment and familiarity differentially shape neural responses to music and highlight the importance of modelling both factors to disentangle their distinct effects on neural activity during music listening.

Deep brain stimulation has demonstrated its therapeutic potential in modulating pathological oscillations associated with Parkinsons disease and epilepsy. However, its efficacy in treating disrupted theta-gamma phase-amplitude coupling seen in memory-related disorders, such as Alzheimers disease, remains poorly understood. While recent studies have targeted the entorhinal-hippocampal circuit, results remain inconsistent. This discrepancy stems from a lack of mechanistic understanding regarding how stimulation protocols affect this circuit. In this work, we present a reduced multicompartment model of the hippocampal CA1 area that reproduces theta-nested gamma oscillations characteristic of healthy neural activity during memory performance. The model comprises pyramidal, basket and OLM cells with simplified morphologies. We also incorporated CA3-to-CA1 axonal projections, providing a foundational framework for studying how stimulation-induced recruitment of afferent pathways modulates CA1 dynamics. By balancing computational efficiency with anatomical accuracy, our model enables systematic investigation of the effects of electrode placement and orientation, as well as stimulation amplitude and frequency on CA1 neural activity. We demonstrate that the excitatory response in CA1 is primarily driven by the recruitment of Schaffer collateral projections. Overall, this work provides a computationally efficient template for exploring diverse stimulation configurations and could be expanded for developing neuromodulatory strategies to restore physiological network dynamics. Author summaryDeep brain stimulation has shown success in treating Parkinsons disease by suppressing abnormal neural activity responsible for movement disorders. However, when applied to memory-related pathologies, such as Alzheimers disease, the therapeutic outcomes remain unpredictable, ranging from cognitive improvement to impairment. This discrepancy highlights a critical gap in our understanding of how stimulation protocols interact with neural dynamics of the targeted circuits. To address this, we developed a computationally efficient model of the hippocampus, which is involved in memory processes, in order to understand how deep brain stimulation might influence its activity. Our model maintains enough biological accuracy to capture essential memory-related neural activity while remaining lightweight enough for rapid execution and systematic exploration of different protocols. This computational efficiency allowed us to conduct systematic investigations of several stimulation configurations to study their effects on hippocampal dynamics. Overall, this model could provide a useful and computationally cost-efficient tool for exploring the mechanisms of deep brain stimulation and help optimize stimulation protocols aimed at alleviating memory disorders.

Individual brains are unique in structure and function. Functional differences are captured by neural fingerprints, which reflect individual differences in behavior and cognition as well as group-level changes related to neurodegenerative diseases. Most research efforts so far have focused on fingerprints com-prising full functional connectomes. However, the high dimensionality of the connectomes can increase computational load and impede performance of machine learning methods in potential applications. A low-dimensional alternative that retains individual features of the full connectomes would thus be beneficial. The present study employed latent-noise Bayesian Reduced Rank Regression (lnBRRR) to learn low-dimensional latent spaces that capture individual features in functional connectivity and power spectral density data derived from MEG recordings. LnBRRR performance was assessed with low training set sizes (N=20-44), and against principal component analysis and linear discriminant analysis. Model performance was also assessed with task data, and the solutions were compared across task conditions with cosine similarity to establish whether individual features are altered by different cognitive processes. LnBRRR captured generalizable individual patterns already at N=20 but N=30-35 was needed to reach optimal test accuracies and to prevent potential overfitting. The model also achieved comparable performance to the alternative models. Latent fingerprints derived from task data attained comparable performance to resting-state latent fingerprints, and lnBRRR solutions were shown to generalize across conditions. Additionally, the model solutions for power spectral density data were discovered to be notably similar, yet differently rotated, over task conditions, suggesting that similar patterns of individual features were captured by the model regardless of the task condition. Altogether, the present results highlight lnBRRR as a potential tool for neuroimaging data analysis and demonstrate that individual differences in power spectral density are largely intrinsic and unaffected by varying cognitive processes.

Intrinsic neural timescales represent the characteristic duration over which information is maintained in neuronal circuits. Evidence suggests that neural timescales vary systematically across the cortical hierarchy, with shorter timescales in primary sensory areas and longer timescales in higher-order association regions. In previous studies, hierarchy has been defined categorically, anatomically, or from the principal gradient of resting-state fMRI functional connectivity derived using diffusion map embedding (DME). Here, we assign hierarchical position to individual human intracranial electroencephalography (iEEG) recording sites by projecting their MNI coordinates onto this embedding space, derived from Human Connectome Project resting-state fMRI data. We estimated neural timescales from resting-state iEEG recordings in adult neurosurgical patients (n=46, 25 female) by extracting the aperiodic component of the local field potential power spectrum using spectral parameterization. Timescales increased monotonically with hierarchical position and associated with two region of interest (ROI)-level measures of network topology derived from DME of participants' iEEG functional connectivity: ROIs with stronger mean functional connectivity exhibited longer timescales, as did ROIs functioning as hubs, defined by proximity to the center of embedding space. Finally, timescales varied with sleep stage, with slowest values during NREM and fastest during wake and REM. The hierarchical gradient present during wake and N1 was no longer detected during REM, N2, and N3 sleep, driven by a selective increase in timescales at lower levels of the hierarchy. This work presents a novel metric of hierarchy that can be applied to iEEG data, establishes a direct link between neural timescales, cortical hierarchy, and network topology in human iEEG, and demonstrates that this hierarchical organization is dynamically modulated by brain state.

In patients with epilepsy, seizures are associated with pathological neural synchronization. However, the preictal period preceding a seizure often exhibits reduced spatial synchronization compared to normal cognition. This observation aligns with the concept of the brain as a complex dynamical system, where a reduction in dimensionality and resilience can precede a phase transition. The Critical Brain Hypothesis suggests a connection between the loss of healthy scale-free behavior and various disorders, including epilepsy. Our study investigates preictal changes by utilizing network features, such as mean node degree and mean clustering coefficient, derived from thresholded correlation matrices of patient intracranial electrocorticographic electrode data. We observed a suppression of intermittent high-synchronization periods within the feature space during the minutes leading up to seizure onset. This constriction of the explored hypervolume in the preictal state indicates a breakdown in the brains ability to maintain normal coherence. We use these preictal changes to predict the probability of seizure onset using a Support Vector Machine algorithm. These discrete predictions can then be combined into real-time continuous seizure risk forecasts via Bayesian updating. This innovative and computationally lightweight approach has the potential to significantly improve upon static predictions, providing opportunities for more adaptable, quantitative, and interpretable tools for managing seizures.

Objective: Aftereffects of 1 Hz repetitive transcranial magnetic stimulation (rTMS) often differ within and between subjects and thus show low reliability. In this study we investigated the mean and individual aftereffects of 1 Hz rTMS using two opposing current directions, their reliability and potential influences of current direction, participants' sex and state on cortical excitability modulations. Methods: Thirteen healthy, right-handed participants underwent four experimental sessions separated by at least 7 days receiving 2000 pulses of suprathreshold 1 Hz rTMS over the primary motor cortex per session. Two sessions were conducted with an induced current direction of anterior-posterior - posterior-anterior (AP-PA) and two sessions with a PA-AP current direction. Before and after rTMS, 100 single TMS pulses were administered with the respective current direction and electromyography was recorded from the first dorsal interosseous. Questionnaires on demographic data and subjective ratings were completed during the experiment. Results: Linear mixed effect model analysis revealed that 1 Hz rTMS induced an excitatory aftereffect when applied with the PA-AP current direction, and no aftereffect with AP-PA. There was a substantial interindividual variability with only three subjects showing an inhibition to 1 Hz rTMS overall. Also, current direction was the only predictor of rTMS aftereffect. Reliability values of these aftereffects were in the poor to moderate range. Conclusions: Current direction plays a crucial role in determining 1 Hz rTMS aftereffects. Reliability was found to be moderate at best. Additional to current direction, more factors need to be considered to tailor the 1 Hz rTMS aftereffects individually.

Brain-computer interfaces (BCIs) have achieved impressive performance by decoding motor and articulatory signals associated with speech production. However, considerably less is known about whether higher-level semantic representations can be decoded from human cortical activity. Demonstrating semantic decoding would advance both our understanding of language organization and the development of BCIs that rely on conceptual rather than purely articulatory information. We recorded intracranial neural activity from patients undergoing stereotactic electroencephalography (sEEG) for clinical epilepsy monitoring while they performed language tasks requiring semantic processing. High-gamma power was extracted from local field potentials and used to generate trial-level features for supervised machine-learning classification. Classification performance was evaluated using cross-validation. Semantic category information was decoded significantly above chance, with mean classification accuracy reaching 29.8% across 15 semantic categories (chance = 6.7%). These findings demonstrate that high-gamma activity contains information about conceptual category membership that can be extracted on individual trials. These results provide evidence that semantic information is accessible from intracranial population recordings and support the feasibility of semantic decoding as a complementary direction for future language BCIs. Beyond neuroprosthetic applications, this work contributes to understanding how conceptual knowledge is represented in the distributed human language network.

Speech cortical mapping by means of navigated repetitive transcranial magnetic stimulation (SCM nrTMS) provides neurosurgeons with noninvasive prior information about individuals cortical speech network. Individualized mapping is required, since the exact locations and activation patterns of speech production show high variability between individuals. We hypothesized that magnetoencephalography (MEG) data of an individuals speech production could guide the SCM TMS process temporally and spatially, leading to higher error rates at MEG-defined locations with TMS pulse timings coinciding with MEG activity. 13 healthy subjects participated in MEG and TMS measurements, where the timing of the TMS pulse (PTI; picture-to-TMS interval) was adjusted based on the individuals MEG activation in a picture naming task. At the group level, significant correlations were observed between the latency of the peak MEG activation and the PTI that produced the highest speech error rate. The MEG peak preceded the best PTI by 132 ms (R=0.713, p=0.006) across the entire stimulation area in the lateral left hemisphere, and by 103 ms (R=0.673, p=0.012) in the left frontal regions. We found 17 combinations of PTI and stimulation area in which the subjects speech error rate increased significantly compared to their average error rate. Our findings suggest that optimal PTIs are highly individual, and that individualizing the PTI according to MEG activation provides a straightforward method for accounting individual variability in speech function and may increase the sensitivity and utility of SCM TMS.

Classification of motor imagery (MI) tasks through EEG is valuable in brain-computer interfacing and rehabilitation engineering. EEG channels selection for MI task classification is well discussed problem and is challenging due to its combinatorial nature. Most of the existing methods are subject and task-dependent. This paper introduces a subject-independent EEG channel selection. The proposed approach consists of two stages. First, we rank channels based on their divergence from a reference channel Cz. We hypothesize that channels less divergent from Cz are more relevant for MI task classification. In the second stage, we employ a three-stage feature selection and classification model to evaluate the selected channels. It consists of a bandpass filter, followed by common spatial pattern (CSP) filter and three classifiers viz. SVM, 1-NN and 5-NN. Two publicly available datasets viz. PhysioNet and BCI Competition III IVa datasets have been used to assess the method. It performs 15.21\% more than 3Cs and just 2.91\% less than all-channels accuracy with as few as 20/118 channels on BCI Competition data and 19.64\% more than 3Cs on the PhysioNet dataset with 16/64 channels. Empirical comparison implies that the method performs better than classical models such as CSP Rank, fishers rank, and normalized mutual information, significantly. Results support that our hypothesis that divergence between channels and a reference channel Cz can be used as a ranking measure for channel selection.

Across contemplative traditions, deeper states of meditation are described as states of heightened clarity, vividness, and stillness of mind, yet what this clarity corresponds to in the brain has remained difficult to specify. The functional signal-to-noise ratio (f-SNR) framework frames mental clarity as a measurable property of neural signals: the degree to which brain activity tracks the causes of sensory signals rather than endogenous, irrelevant fluctuations. It predicts that deepening meditation should raise f-SNR, expressing sensory events more faithfully in neural signals against ongoing background activity. We tested this prediction across different levels of meditative depth. Twenty-nine experienced Vipassana practitioners meditated while auditory tones were presented, periodically reporting their depth of meditation. f-SNR was quantified from event-related potentials (ERPs) in a fronto-central P3 window and from single-trial decodability of auditory tone-evoked activity against no-tone background EEG. High-depth states were associated with greater ERP signal-to-noise ratio, stronger single-trial signal consistency, and improved decodability of auditory tones. These results suggest that meditative depth is expressed in the reproducibility and stimulus-background separability of sensory responses, consistent with deep meditation enhancing the brain's functional signal-to-noise ratio by improving the clarity of sensory signals and reducing endogenous noise.

The beta oscillations of the human primary sensorimotor cortex (SM1) play a crucial role in regulating motor and cognitive behavior in health and disease. However, their assessment relies on costly and complex neuroimaging techniques, limiting scalability and translational applications. We present a novel method for assessing beta oscillations from easily obtained peripheral electromyography and force recordings. We show that movement-induced modulations in SM1 beta oscillations can be assessed from the electromyography or force recordings of a contracted contralateral hand muscle. We demonstrate the fidelity of this method in young and elderly healthy participants and in Parkinsons disease patients. We also demonstrate that the resting-state SM1 beta interhemispheric coupling can be assessed from an interhand coupling between the electromyography or force of contracted homologous hand muscles. This methodology enables scalable and cost-effective investigations of beta oscillations for all fields of human neuroscience and for the development of accessible disease/therapeutic markers.